o 
o 



> 
oo 
o 
o 
m 
o 

o 

"a 

• I— I 

> 
X 



Exact soliton solutions, shape changing collisions and partially coherent solitons in 

coupled nonlinear Schrodinger equations 
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We present the exact bright one-sohton and two-sohton solutions of the integrable three coupled 
nonlinear Schrodinger equations (3-CNLS) by using the Hirota method, and then obtain them for the 
general A^'-coupled nonlinear Schrodinger equations (A'^-CNLS). It is pointed out that the underlying 
solitons undergo inelastic (shape changing) collisions due to intensity redistribution among the 
modes. We also analyse the various possibilities and conditions for such collisions to occur. Further, 
we report the significant fact that the various partial coherent solitons (PCS) discussed in the 
literature are special cases of the higher order bright soliton solutions of the A^-CNLS equations. 

PACS numbers: 42. 81. Dp, 42.65. Tg 



In recent years the concept of soliton has been receiving 
considerable attention in optical communications since 
soliton is capable of propagating over long distances with- 
out change of shape and with negligible attenuation [1-3] . 
It has been found that soliton propagation through opti- 
cal fiber arrays is governed by a set of equations related 
to the CNLS equations[l,2]. 



N 

p=i 



iqjz+qjtt + 2fi\^\qp\^qj ^0, J = 1, 2, ...TV, (1) 



where qj is the envelope in the jth core, z and t rep- 
resent the normalized distance along the fiber and the 
retarded time, respectively. Here 2fi gives the strength 
of the nonlinearity. Eq. (1) reduces to the standard enve- 
lope soliton possessing integrable nonlinear Schrodinger 
equation for N = 1. For N = 2, the above Eq. (1) 
governs the integrable Manakov system [4] and recently 
for this case the exact two-soliton solution has been ob- 
tained and novel shape changing inelastic collision prop- 
erty has been brought out [5]. However, the results are 
scarce for A'^ > 3, even though the underlying systems 
are of considerable physical interest. For example, in 
addition to optical communication, in the context of bio- 
physics the case = 3 can be used to study the launch- 
ing and propagation of solitons along the three spines of 
an alpha-helix in protein [6]. Similarly the CNLS Eq. (1) 
and its generalizations for iV>3 are of physical relevance 
in the theory of soliton wavelength division multiplex- 
ing[7], multi-channel bit-parallel- wavelength optical fiber 
networks[8] and so on. In particular, for arbitrary N, 
Eq. (1) governs the propagation of A^-self trapped mutu- 
ally incoherent wavepackets in Kerr-like photorefractive 
media[9] in which qj is the jth component of the beam, 
z and t represents the normalized coordinates along the 
direction of propagation and the transverse coordinate, 
respectively, and kpP represents the change in the 

refractive index profile created by all incoherent compo- 
nents of the light beam [9] when the medium response is 



slow. The parameter n = fcQn2/2, where n2 is the nonlin- 
ear Kerr coefficient and kg is the free space wave vector. 

In this letter, we report the exact bright one and two 
soliton solutions, first for the N — 3 case and then for 
the arbitrary N case, where the procedure can be ex- 
tended in principle to higher order soliton solutions, us- 
ing the Hirota bilinearization method. In particular, we 
point out that the shape changing inelastic collision prop- 
erty persists for the > 3 cases also as in the N — 2 
(Manakov) case reported recently [5] , giving rise to many 
possibilities of energy exchange. Furthermore, we point 
out that in the context of spatial solitons the partially 
coherent stationary solitons(PCS) reported in the recent 
literature[9-10] are special cases of the above general soli- 
ton solutions which undergo shape changing collisions. 

The bright one-soliton and two-soliton solutions of the 
3-CNLS system, 

iqj,+qjtt+2^i{\qi\^ + \q2\^ + \q3\^)q,=0, J = 1,2,3, (2) 

can be obtained from its equivalent bilinear form result- 
ing from the transformation qj = g'^-* //, 



0, D^f.f^2^i 



3 

n=l 



(3) 



where * denotes the complex conjugate, Dz and Dt are 
Hirota's bilinear operators [11], and g^^^^s are complex 
functions, while f{z,t) is a real function. The resulting 
set of Eqs. (3) can be solved recursively by making the 
power series expansion g^^^ = X5i"'^ +X^.93"'^ + ■ • / = 1 + 
X^/2+X^/4+- ■ -J J = 1: 2, 3, where x is a formal expansion 
parameter. In order to get the one-soliton solution, the 
power series expansions are terminated as g^^^ = ^'"^■'^ 



and / = 1 + After deducing t;'-') and f as g 



xgi 

U) = 



1,2,3 and / = 1 



gVi+Vi+R where 



fi J2^j=i I'^i"'' 1^/ i^i + ^1)^1 the bright one-soliton solution 



is obtained as 



(91,92,93) 



/ (1) (2) (3)V 

\^a\ ',al ',a\ ' j 



1 



cosh [r]iR + f j 



(4) 



where r]i=ki(t + ikiz), Aj=aY / A and A = 

(m(E-=i Here a['\ fci, j = 1,2,3, are four 

arbitrary complex parameters. Further kmAj gives the 
amplitude of the jth mode and 2fci/ the soliton velocity. 

The general bright two-soliton solution of Eq. (2) 
can be generated by terminating the series as g^^^ = 
X9i^ + X^ga^ and / = 1 + x'^f2 + x'^fi and solving the 
resultant linear partial differential equations. This so- 
lution contains eight arbitrary complex parameters, a'p 
and ki, I = 1,2 and j = 1, 2, 3. It is given by 



(5) 



where £) = 1 + e^^+'^'+^^ + e''i+''2+'5o + Qnl+m+^o + 



-Ri 



Kii/(fci + kl), e"-' = K22/(fc2 + k^), = «;i2/(fci + k^), 



-k2){a[^^K2i- 
.,0) 



,.0') 
2 

0), 



Kii))/((fci+A;i*)(A;i* + fc2)), 



e*2.- = ((fc2-fcl)(a^^^M2-<V22))/((fc2+fc^)(fcl+fc|)), 

(|fcl-fc2p(KllK22-Kl2«;2l))/((fcl+fci*)(fc2+fc5)|fci + 



0^3 



fc2*n and Kii = I^EUiC^i4 /iki + kt), i, 1 = 1,2, 
j = 1,2,3. Though one can proceed to obtain higher 
order soliton solutions in principle by making use of the 
general power series expansion, the details become com- 
plicated and we will present the results separately. 

The nature of the interaction of the underlying soli- 
tons can be well understood by making an asymptotic 
analysis of the two-soliton solution [5]. Asymptotically, 
the two-soliton solution(5) can be written as a combina- 
tion of two one-soliton solutions and their forms in the 
two different regimes 2; ^ — 00 and z ^ 00 are similar 
to those of the one-soliton solution given in Eq. (4) but 
differing in amplitude (intensity) and phase. The analy- 
sis reveals that there is an intensity exchange among the 
three components of each soliton during this two-soliton 
interaction, which can be quantified by defining a tran- 
sition matrix Tj such that A^+=A^~T^-, j=l,2,3 and 
Z=l,2, where the superscripts l± represent the solitons 
designated as and ^2 at 2; — > ±00 , and denote 
the corresponding amplitudes. 

Consequently, the three modes qi , 52 and (73 of Si hav- 
ing magnitudes of amplitudes |A^~|fci_R ^\a^i^\kifj/ Ai, 
where Ai={ii{J2^^-^ \a'"/^ P))^^^; exchange intensity given 
by the square of the transition matrices, \Tj\'^ = |1 — 

A2(a2''Vai''^)l^/|l - ^1^2!, j = 1, 2. 3, along with a phase 
shift $^ = (i?3 - i?2 - -Ri)/2 during collision. Here 
Ai = K21/K11 and A2 = ki2/k22. In a similar fashion 
due to collision the three modes qi , (72 and q^ of 5*2 also 
exchange an amount of intensity, = |1 — A1A2I/II — 

Xi{a'"/^ /a2 ^)\'^, j = 1,2,3, respectively and change their 
amplitudes to \A'^.~^\k2R = l^a''^ |fc2K/A2 from \A^f\k2R, 



respectively Here A2=(/i(5];^^ Ja^-'^ P))^/^ The asso- 
ciated phase shift for this soliton is = — (i?3 — R2 — 
R\)/2. We also note there is a net change in the relative 
separation distance between the solitons due to collision 
by Aa;i2 = (/cia + k2R)\{Rz - R2 - Ri)\/2kiRk2R. 
Also, we note that for the special case \T- \ = 1, 1 = 1,2, 

j = 1,2,3, which is possible only when a^^/a^^ = 
a^^^/a^2^= af^/a'"2 \ ttic collision corresponds to the 
standard elastic collision. For all other cases, the quan- 
tity \Tj\ ^ 1, which corresponds to a change in the val- 
ues of the amplitudes of the individual modes leading 
to a redistribution of the intensities among them and 
corresponding to a change in the shape of the soliton. 
However, during the interaction the total intensity of the 
individual solitons 5*1 and 5*2 remains conserved, that is 

|4^|2 + |4T|2 + |Af |2 = l/^. 

The above shape changing (inelastic) collision during 
the two-soliton interaction of the 3-CNLS can occur in 
two different ways. The first case is an enhancement 
of intensity in anyone of the modes of either one of the 
solitons (say ^i) and suppression in the remaining two 
modes of the corresponding soliton with commensurate 
changes in the other soliton 82- The other possibility is 
an interaction which allows one of the modes of either 
one of the solitons (say Si) to get suppressed while the 
other two modes of the corresponding soliton to get en- 
hanced (with corresponding changes in ^2). In either of 
the cases, the intensity may be completely or partially 
suppressed(cnhanced). Thus as a whole during the in- 
elastic interaction among the two one solitons and 
52 of the 3-CNLS, the soliton Si {S2) has the follow- 
ing six possible combinations to exchange the intensity 
among its modes: {qi,q2,q3) ^ (g?, 92- 93)11 («> c = 
5 (suppression), (enhancement)) with i = 1, a = E, 
b = S, c = S; i = 2, a = S, b = E, c = S; i = 3, a = S, 
b = S, c = E; i = 4, a = S, b = E, c = E; i = 5, a = E, 
b = S, c = E andi = 6, a = E,b = E, c = S. 

Two of the above interactions involving a dramatic 
switching in the intensity are depicted in Fig. 1 for illus- 
trative purpose for specific choice of soliton parameters. 
These may also be viewed as the two-soliton interaction 
in a waveguide supporting propagation of three nonlin- 
ear waves simultaneously. For other choices, in general, 
partial suppression(enhancement) of intensity among the 
components will occur depending on the values of the 
transition matrix elements Tj. Fig. la is plotted for the 

parameters ki = I + i, k2 = 2 — i, a'^^ = a^^ = 
(39 + i80)/89, a'-^^ = af^ = af^ = = 1 and 

/i = 1. In this figure the intensities of the components qi 
and q2 of 5*1 (^2) are almost completely suppressed (en- 
hanced) and that of the third component is enhanced 
(suppressed). The second possibility of enhancement 
(suppression) of intensity in the qi and 52 components of 
Si{S2) and suppression(enhancement) of intensity in its 



2 



(73 component are shown in Fig. lb, in which the parame- 
ters are chosen as fci = l + i, ^2 = 2 — i, a^^^ = 0.02 + O.li, 
(2) „ 1 • (3) 



X9^^ + X'g^' + X^g^' and / = 1 + x'f2 + XV4 + xVe, 



-.(2) 



.(3) - 



= 1. 



Now it is straight forward to extend the above analysis 
to obtain the one-soUton and two-sohton solutions of the 
arbitrary A'^-CNLS Eq. (1). After making the bilinear 
transformation qj = g^^^ //, j = 1, 2, 3, . . . , A'' in Eq. (1), 
one can get a set of bilinear equations of the form (3) 
but now with j, n = 1, 2, 3, . . . , A. Then by expanding 
(/(■'^s and / in power series up to A^terms and following 
the procedure mentioned above, the one-soliton and two- 
soliton solutions of Eq. (1) can be obtained. 

(a) One-soliton solution: 



{qi,q2,...,qN) 



kiRe 



cosh(?7ifl + f ) 



{Ai,A2,...,ANy , (6) 



where r/i = ki{t + ikiz), Aj = af^A, A = 

(M(Ef=i e« = AV(fci + a^P and fci, 

j=\, 2 . . . , N, arc (A^ + 1) arbitrary complex parameters. 

(b) Two-soliton solution: This solution will also be 
of the same form as Eq.(5) with the replacements, j = 



1,2, 



, A^ and k 



ii 



(n) (n)* 



a 



/{h 



k^), where 



i, I = 1,2. One can also verify that this two-soliton solu- 
tion will depend on 2(A^-|- 1) complex parameters and the 
shape changing interaction can lead to intensity redistri- 
bution among the modes of the soliton of the A^-CNLS 
system in 2-'^ — 2 ways (by generalizing the N = 3 case). 
We believe that such studies will have important appli- 
cations in logic gates and all optical computations[12]. 

From an application point of view, it has been ob- 
served recently that the CNLS equations can support a 
kind of stationary solutions known as partially coherent 
solitons (PCS). In particular, explicit form of such solu- 
tions have been given for A^ = 2,3 and 4 cases of Eq. 
(1) in Ref.[9]. They have also been shown to have vari- 
able shapes. Now, the generalized Manakov equation(l) 
is intcgrable[13] and hence its 7V-soliton solution can be 
obtained in principle by extending our above analysis. 
So the natural question arises as to what is the relation 
between the PCS and the exact A'-soliton solutions. 

To answer the above question, let us look at the 
A' = 2,3 and 4, cases of Eq. (1) explicitly. One can 
check that very special cases corresponding to specific 
parametric restrictions in the two-soliton solution of the 
N = 2 case, the three-soliton solution of the N = 3 case 
and four-soliton solution of the A = 4 case give rise to 
the 2-soliton, 3-soliton and the 4-soliton PCSs, respec- 
tively, reported in [9] . In order to appreciate this we con- 
sider as an illustration the three-soliton solution of the 
N = 3 case of Eq. (1). Instead of writing down the full 
3-soliton solution of the N = 3 case explicitly and choos- 
ing the special parametric values, we make the follow- 
ing simplified procedure. Starting from the bilinear Eqs. 
(3) and terminating the series for g^^^ and / as g^^^ = 



(i) (i) 
one can identify gl = a\ e^^ 



Si), 



e^^ where 



Vn = kn{t + iknz), j, n = 1, 2, 3 in which al and ki are 
complex parameters. Finding g^\ g^\ j = 1, 2, 3, /2, 
and /e, the three soliton solution is obtained. Instead, 
as a special case, we look for a stationary solution with 
= 0, a« = = ) = af = af^ = = 

and o'l^ = —(^^2^ = Oi'i'^ = 1, in order to gain insight 
into the physics of the problem. Then, the resulting ex- 
plicit expression for the three-soliton solution has been 
found after simple algebraic manipulation to be exactly 
the same as the stationary PCS for N=3 given in Eq.(17) 
of Ref.[9a]. One can also check that with the choice 

kni = 0, 02^ = a^i^ = 0, a^^-* = —a^'^ = 1, in the 
two-soliton solution of the N = 2 case (Manakov equa- 
tion) [5] of Eq. (1), the Af = 2 soliton complex(PCS) is 
obtained. By a similar analysis we have verified that the 
A" = 4 PCS also results as a special case of the 4-soliton 
solution of the 4-CNLS equation. Extending this idea it 
is clear that the PCS which is formed due to a nonlinear 
superposition of A"- fundamental solitons [9] is a special 
case of Af-soliton solution of the A^-CNLS Eq. (1). 

Further, it has been found that these PCS are of vari- 
able shape and also change their shape during collision 
with another PCS [9]. The reason for the shape change of 
the PCS can be deduced from the interaction properties 
of the solitons discussed above. The solitons are char- 
acterized by their amplitudes AjkiR and their velocities 
2fc;7(so that the angle of incidence is Oi = tan~^ {2kii)) . 
During a pair-wise interaction of two fundamental soli- 
tons of A^-CNLS equation there is an energy sharing be- 
tween them resulting in a novel shape changing colli- 
sion, depending on the transition matrix elements Tj, 
the phase shift and change in the relative separa- 
tion distance Axij defined earlier . Since the PCS is 
a special case of the A^-soliton solution, parametrized as 
above, it naturally possesses a variable shape. For exam- 
ple in Figs.l, the solitons Si and 52 are travelling with 
velocities 2kij = 2 and 2^2/ = —2, respectively. For 
the chosen parameters, in Fig. la, the amplitudes of the 
modes of the solitons Si and ^2 before interaction given 
respectively by kmlA]'] = (0.577,0.577,0.577) and 



Aj 



(0.857,0.857,1.591) change to kiii]AY 



(0.183,0.183,0.966) and fc2fi|A|+| = (1.155,1.155,1.155) 
after interaction, preserving the total intensity of 
each of the soliton. Similarly in Fig. lb, the ampli- 
tudes of the solitons Si and 5*2 before interaction are 
(0.101, 0.099, 0.990) and (1.322, 1.335, 0.686) respectively, 
while after interaction they become (0.466,0.484,0.741) 
and (1.155, 1.155, 1.155). The phase shifts suffered by the 
solitons are $1 = -^^ = -0.787(Fig.la),-0.600(Fig.lb). 
During the collision process the initial separation dis- 
tance — —0.668 changes to = 0.513 in Fig. la and 
-0.036 to 0.865 in Fig.lb. 
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(a) (b) 
FIG. 1. Intensity profiles of tlie tliree modes of tlie 
two-soliton solution in a waveguide described by the CNLS 
Eq. (2) showing two different dramatic scenarios of the shape 
changing collision. The parameters are chosen as in the text. 

In a similar way, the variable shape of the PCS during 

interaction with another PCS also arises from the fun- 
damental bright soliton collision of the Manakov system. 
The collision of two PCS each comprising of m and n 
soliton complexes respectively, such that m + n = A'^, 
is equivalent to the interaction of N fundamental bright 
solitons (for suitable choice of parameters) represented 
by the special case of A^-soliton solution of the A^-CNLS 
system. Further details will be published elsewhere. 

Our above analysis has considerable practical relevance 
in view of the various recently reported interesting ex- 
perimental observations. Firstly, the Manakov spatial 
solitons have been observed in AlGaAs planar waveg- 
uidcs[14] and their collisions involving energy exchange 
of precisely the type discussed here has been experimen- 
tally demonstrated[15]. Also collisions between PCS's of 
shape changing type as treated here were observed in a 
photorefractive strontium barium niobate crystal using 
screening solitons[16]. Further partially incoherent soli- 
tons have been observed through excitation by partially 
coherent light[17] and with a light bulb[18]. Using dif- 
ferent techniques, such as the coherent density function 
theory [19], to describe such incoherent solitons one can 
obtain the N-coupled NLS equations of the form(l) con- 
sidered in this Letter. We believe that our exact analyt- 
ical results will give further impetus in the experimental 



investigations of these solitons. 

In conclusion, we have shown that iV-CNLS Eq. (1) 
possesses fascinating type of soliton solutions undergo- 
ing shape changing (inelastic) collision property due to 
intensity redistribution among its modes. The many pos- 
sibilities for such collisions to occur provides interesting 
avenues of research in developing logic gates and in all- 
optical digital computations[12]. We have also shown the 
interesting fact that the multisoliton complexes are spe- 
cial cases of the shape changing bright soliton solutions 
and pointed out that the variable shape of PCS is due to 
the shape changing collision of the fundamental solitons 
which is an inherent nature of the iV-CNLS system. 

The work of M.L. and T.K. forms part of a Dept. of 
Science and Technology, Govt, of India research project. 
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